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High-redshift QSOs in the FIRST survey 
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ABSTRACT 

In a pilot search for high-redshift radio QSOs, we have obtained spectra of 55 FIRST 
sources (Siaghz > 1 mJy) with very red (0-E > 3) starlike optical identifications. 10 
of the candidates are QSOs with redshifts 3.6 < z < 4.4 (4 previously known), six with 
z > 4. The remaining 45 candidates comprise: a z — 2.6 BAL QSO; 3 low-redshift 
galaxies with narrow emission lines; 18 probable radio galaxies; and 23 M stars (mainly 
misidentifications) . The success rate (high-redshift QSOs / spectroscopically-observed 
candidates) for this search is 1/2 for Siaghz > 10 mJy, and 1/9 for Siaghz > 1 rnJy. 
With an effective search area of 4030 deg 2 , the surface density of high-redshift (z > 
4) QSOs discovered with this technique is 0.0015 deg -2 . 
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1 INTRODUCTION 

Observations of high-redshift (z > 4) QSOs strongly con- 
strain models of: the formation and evolution of galaxies 
and their central black holes (Kauffmann & Haehnelt 2000) ; 
the contribution of QSOs to the ionisation of the intergalac- 
tic medium at high redshifts (Steidel, Pettini, Adelberger 
2001); and, through studies of the Lyman absorption forest, 
the chemical evolution of the intergalactic medium along the 
line of sight to the QSO (Rauch 1998, Hamann & Ferland 
1999). 

Most of the ~ 300 z > 4 QSOs now known have been 
discovered as a result of searches for objects with unusually 
red optical colours (Fig. 1). A summary of recent searches is 
given in Table 1. Kennefick et al (1995a, b) searched for ob- 
jects with unusually red g, r, i colours in 340 deg 2 covered by 
the second Palomar Sky Survey (POSS-II) of the northern 
hemisphere and discovered 10 z > 4 QSOs (see also Djor- 
govski 2001). Similarly, Sloan Digital Sky Survey (SDSS) 
commissioning images have been searched for objects with 
red u,g,r,i,z colours. More than 100 z > 3.5 QSOs have 
been found so far (Fan et al 2000, 2001; Schneider et al 2001; 
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Zheng et al 2001; Anderson et al 2001, Richards et al 2001), 
and this number will rise rapidly as SDSS nears completion 
(~ 10000 deg 2 ) over the next few years. In the south, Storrie- 
Lombardi et al (2001) searched UKST plates scanned with 
the Automated Plate Measuring facility (APM) in Cam- 
bridge for objects with Bj — R > 2.5, and found 49 QSOs 
with z > 4. 

A disadvantage of the purely optical searches is that 
complex colour criteria are needed to exclude the much 
larger number of red stars, which makes it difficult to correct 
for incompleteness. Alternatively, simple one-colour criteria 
can be used, but then a large fraction of the high-redshift 
QSOs is missed (compare in Table 1 the surface densities at- 
tained for different types of optical selection) . Starting with 
a sample of radio QSOs allows one to reduce the number 
of candidates without resorting to complex colour criteria, 
and is less likely to bias the selection against dusty objects. 
E.g. Hook et al (1998) sought red objects identified with 
flat-spectrum QSOs with Ssghz > 25 mJy over 1600 deg 2 , 
and found 6 with z > 3 (none with z > 4). Snellen et al 
(2001) sought red objects identified with flat-spectrum radio 
sources S$ghz > 30 mJy in the Cosmic Lens All-Sky Sur- 
vey (CLASS, Myers et al 2001) over 6400 deg 2 , and found 4 
QSOs with z > 4, i.e. 1 z > 4 QSO per 1600 deg 2 . 

Higher surface densities can be attained by observing 
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Redshift z 

Figure 1. The variation of POSSI O — E colour with redshift, 
for QSOs from the FIRST bright quasar survey (crosses, White et 
al 2000), the SDSS and POSSII searches for high-rcdshift QSOs 
(spots), and a selection of radio QSOs with 3.4 < z < 4.0 from 
the NASA Extragalactic Database (circles). The QSOs in these 
samples were selected in part on the basis of colour, but nev- 
ertheless trace the typical reddening of a QSO with increasing z 
(prediction for typical QSO shown as solid line) as the continuum 
redward of the Lya line moves out of the observed O band (red 
limit 5000 A). As redshift increases above 4.5, the predicted 
brightness of a QSO in E band drops sharply, so few are likely to 
be detected. With the criteria O — E > 3, E < 18.6 used here, 
our QSO search is probably complete for z > 4, but incomplete 
for 3 < z < 4. 



E mag 

Figure 2. Distribution in colour and magnitude of starlikc 
(POSSI/ APM) optical identifications of FIRST radio sources. 
Objects with no O magnitude have been plotted with O — E = 
21.7 - E, and the E magnitudes have been randomised ± 0.1 so 
that the density of points in this region can be seen. Most of the 
objects in the lower part of the diagram are true QSOs. Most of 
those clustered at upper right are misclassified radio galaxies (the 
dashed line shows the expected colour-magnitude relationship for 
a redshifted giant elliptical galaxy), and many of the other red 
objects arc misidentifications with galactic M stars. Some may 
also be BAL QSOs, which tend to be red. The candidate high- 
rcdshift QSOs observed here are drawn from the region enclosed 
by the dotted line. 



the counterparts of fainter radio sources. In this paper we 
report on a pilot search for high-redshift QSOs amongst very 
red, O — E > 3 (see Fig. 1) , optical identifications of FIRST 
radio sources Siaghz > 1 rnjy (with no selection on radio 
spectral index). 

There is no overlap between this search and the FIRST 
bright quasar survey (White et al 2000, Becker et al 2001), 
whose selection criterion O — E < 2 preclude detection of 
most high-redshift QSOs (none with z > 3.7 were cata- 
logued), as do the colour selection criteria of most large QSO 
surveys, e.g. the 2dF QSO survey (Croom et al 2001). 



2 SAMPLE 

The FIRST radio catalogue (White et al 1997) currently in- 
cludes 722354 sources detected with peak Siaghz > 1 mjy 
over 7988 deg 2 (» 90 sources deg~ 2 ), mainly 7 < RA < 17 h , 
-5 < Dec < 57°. We sought optical identifications of these 
sources with objects catalogued by the APM (Irwin et al 
1994) on the POSS-I O (blue) and E (red) plates. We did 
not seek optical identifications at the mid-points of likely 
double radio sources (~ 10% of the catalogue), since this 
would have required a substantial enlargement of the search 
area. The colour-magnitude distribution of the star-like opti- 
cal identifications lying within 1.5 arcsec of a FIRST source 



is shown in Fig. 2. We selected as candidate z > 4 QSOs 
all starlike optical identifications with O — E > 3 (see Fig. 
1) and E < 18.6 (the POSS blue limit is O ~ 21.7). Trial 
spectroscopy showed that ~ 80% of these are actually radio 
galaxies, which greatly outnumber QSOs at E ~ 18, O — E ~ 
3 (Fig. 2). It was possible to filter out most of these by in- 
specting images from the digitised POSS-II survey, using the 
SExtractor program (Bertin & Arnouts 1996) to characterise 
the relationship between FWHM and intensity for stars (in 
each image), and thus to distinguish between stellar and ex- 
tended images. In addition, we filtered out the » 25% of 
starlike images for which APM recorded no blue magnitude 
(i.e. implied magnitude fainter than the blue plate limit), 
but which were easily visible on the Minnesota APS scans 
of POSS-I (Pennington et al 1993); these images were missed 
by APM due to the coarser scanning resolution and greater 
susceptibility to confusion with nearby images. This selec- 
tion procedure yielded a total of 109 candidate high-redshift 
QSOs. 



3 OBSERVATIONS AND REDUCTION 

Spectra were obtained of a random sample of 55 of the 109 
candidates with the Isaac Newton Telescope's IDS spectro- 
graph on 2001 Mar 14, and with the Isaac Newton (IDS), 
William Herschel (ISIS) and Calar Alto 2.2-m (CAFOS) 
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Table 1. Recent searches for high-redshift QSOs 



Search programme 


Colour 


optical 




Area 


Cands 


z > 


N 


<7 


Hiz / 


Reference 




criterion 


limit 


mjy 


deg 2 








deg -2 


cands 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


DPOSS-II 


gri 


r < 19.6 


_ 


340 


85 


4.0 


10 


0.029 


0.11 


Kennefick et al 1995a, b 


SDSS 


ugriz 


i < 20 




180 




4.0 


18 


0.10 




Fan et al 2001 


APM/UKST 


Bj - R >2.5 


R < 19.3 




8000 




4.0 


49 


0.0061 




Storrie-L. et al 2001 


INT wide-angle surv. 


griz 


r < 24 




12.5 




4.0 


3 


0.24 




Sharp et al 2001 


CLASS fiat-spectrum 


O - E > 2 


E <19 


S s > 30 


6400 


27 


4.0 


4 


0.00063 


0.15 


Snellen et al 2001 


GB/FIRST flat-spec. 


O - E > 1.2 


£ <19.5 


S s >25 


1100 


50 


4.0 





0.0 


0.00 


Hook et al 1998 




O - E > 1.2 


E <19.5 


S 5 >25 


1100 


50 


3.0 


6 


0.005 


0.12 


u 


FIRST 


O - £ > 3 


E <18.6 


51.4 >1 


4030 


55 


4.0 


6 


0.0015 


0.11 


This paper 




O - £ > 3 


£ <18.6 


51.4 >10 


4030 


9 


4.0 


4 


0.0010 


0.44 





Column 4 gives, where appropriate, the radio flux-density limit S v , for frequency v in GHz. Column 9 gives the surface density 
of high-redshift QSOs discovered. Column 10 gives the ratio between the number of high-redshift QSOs found (column 8) and the 
number of candidates for spectroscopy (column 6). 




E mag 

Figure 3. Colour-magnitude distribution for the candidate high- 
redshift QSOs. The symbols for spectral type are as used in Ta- 
ble 2 (Q = QSO, S = narrow-emission-line galaxy, G = radio 
galaxy, * = star, ? = probable radio galaxy). Large font indicates 
Si.4GHz > 10 mjy; small font Si^GHz < 10 mjy. A nominal 
POSSI blue plate limit O = 21.7 is indicated by the dashed line. 
Objects with no O magnitude have been plotted with O — E = 
21.7 - E, and the E magnitudes have been randomised ± 0.2 mag 
to reduce crowding. The arrows indicate upper or lower limits. 



spectrographs on service nights between 2000 March 19 and 
2001 Jan 5 (dates given in Table 2). Apart from 2000 Apr 5, 
May 12 and 2001 Jan 5, the nights were photometric. The 
INT IDS spectrograph was used with the R150V grating, 
yielding spectra with dispersion 6.5 A/pixel, usually cen- 
tred at 6500 A. The WHT ISIS dual-arm spectrograph was 
used with the R158R and R158B gratings in the red and 
blue arms respectively (with a dichroic separating the light 
blue and red of 6100 A). In the red arm, with a TEK CCD 
detector, this yielded spectra with dispersion 2.9 A, usually 
centred at 7400 A. In the blue arm, with an EEV CCD, the 
dispersion was 1.6 A/pixel, and the spectra were usually 



centred at 4800 A. The CAFOS spectrum has dispersion 8.5 
A /pixel. 

The images were debiased and flat-fielded, cosmic 
rays were eliminated, and the spectra were extracted, 
wavelength-calibrated and intensity-calibrated in the usual 
way, using the IRAF package. A summary of the observa- 
tions and results is given in Table 2. 



4 RESULTS 

Spectra showing clear emission or absorption features were 
classified 'Q' (QSO, broad emission lines), 'S' (Sy2 or star- 
burst type spectrum, narrow emission lines) or '*' (spectral 
features of M star). The remainder have been classified 'G' 
(probable radio galaxy, i.e. giant elliptical), where the char- 
acteristic 4000-A break is detected, otherwise '?'. Most of 
the '?' objects will be radio galaxies, rather than misiden- 
tifications with stars, judging from the number of M stars 
detected, the typical distribution of spectral types at this 
magnitude, and the magnitude distributions and the distri- 
butions on the sky of the '?' and '*' objects. One of the '?' 
objects, 1036+42, has a flat SED inconsistent with it be- 
ing a radio galaxy; it may be a low-redshift QSO (z < 2.2). 
Of the other '?' objects, none has an SED consistent with 
it being a high-redshift QSO, except perhaps 1120+34 and 
1349+38. The redshifts of the high-redshift QSOs were mea- 
sured from emission lines other than Lya (because of the 
asymmetrical absorption), except in the case of 0839+51. 
The colour-magnitude diagram of the candidates is shown 
in Fig. 3. The spectra of the high-redshift QSOs are pre- 
sented in Figs. 4, 5. 

We find 10 high-redshift QSOs (z > 3.6) out of 55 candi- 
dates. These candidates were selected from 109 for the whole 
FIRST catalogue, covering 7988 deg 2 , so the effective area 
covered by this search is ~ 4030 deg 2 , i.e. we find 0.0025 
high-redshift QSOs deg" 2 with z > 3.6, or 0.0015 deg" 2 
with z > 4.0. For Siaghz > 1 mjy, E < 18.6, the z > 4 
search will be near-complete (Fig. 1). The surface density 
on the sky is compared with that from other searches for 
high-redshift QSOs in Table 1. The surface density of z > 4 
radio QSOs discovered here is twice as high as that from the 
Snellen et al search, because of the lower flux-density limit. 
The efficiency of the search (number of high-redshift QSOs / 
number of candidates for spectroscopy) is 0.44 for Siaghz > 
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Wavelength (A) Wavelength (A) 

Figure 4. Spectra of 8 of the 10 high-rcdshift QSOs (sec Fig 5 for the remainder). Spectral features are labelled at wavelengths 
corresponding to the quoted redshift, assuming rest-frame wavelengths in A of 1216 (Ly«), 1240 (NV), 1302 (Ol/Sill blend), 1400 
(SilV/OIV] blend), 1549 (CIV). For 0941+51, the metal-line absorption systems at z = 3.80, 3.83 are also indicated. The spectra have 
not been corrected for the terrestrial atmospheric absorption bands, notably at 7594 and 6867 A (A and B bands). 
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Table 2. High-rcdshift QSO candidates 



RA 


Dec 


Si. 4 


R-O 


E 


O-E 


Tel 


Date 


T 


ID 


z 




Notes 


J2000 


J2000 


mjy 


ii 








obs. 


sec 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


High-redshift QSOs: 
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524517.1 
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3 Q 
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3 87 
O.o / 


03 

.uo 
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083946.11 
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41.6 


0.8 


18.5 


>3.2 


J 


010314 


300 


o 


4.42 


.02 


Snellen ct al (2000), z = 4.41 


094119.36 


511932.3 


2.7 


1.3 


17.6 


>4.1 


w 


000319 


1200 


Q 


3.85 


.01 




105756.27 


455553.0 


1.4 


0.5 


16.5 


3.9 


J 


000405 


300 


Q 


4.11 


.02 


Stern et al (2000), z = 4.12 


122027.93 


261903.6 


35.0 


0.4 


17.9 


3.3 


J 


010314 


1200 


Q 


3.694 


.005 




130940.61 


573309.1 


11.5 


1.1 


18.1 


>3.6 


w 


000319 


300 


Q 


4.274 


.005 




142308.19 


224157.5 


35.4 


0.9 


18.5 


>3.2 




010314 


1200 


Q 


4.316 


.005 




151601.55 


430931.3 


1.5 


0.9 


18.4 


>3.3 


I 


010314 


300 


Q 


2.590 


.005 




163950.51 
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23.8 


0.5 


17.1 


3.9 


w 


000718 


1000 


Q 


3.99 


.015 




Low-redshift QSOs and galaxies: 
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The columns give: (1-2) right ascension and declination, (3) FIRST 1.4-GHz (integrated, sometimes < peak) flux density, (4) radio 
- optical displacement, (5-6) APM POSSI E magnitude and O — E colour, (7) telescope with which spectrum obtained (W = 
William Herschel 4.2-m, I = Isaac Newton 2.5-m, C = Calar Alto 2.2-m), (8) date of observation, (9) exposure time, (10) spectral 
type (Section 4), (11-12) redshift and rms error on redshift. For low-redshift emission-line objects, the rest-frame wavelengths of 
detected lines arc given in A in the last column: 3727 Oil, 4861 H/3, 4959/5007 OIII, 6563 Ha. For high-redshift QSOs, see Figs. 
4, 5. 
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Figure 5. Spectra of the two remaining high-rcdshift QSOs (see caption of Fig. 4); and of 4 other emission-line objects detected in this 
search (see text). 1516+43 was undetected on the POSS blue plate, but the field is optically confused and its true colour may be bluer 
than O - E = 3. 



10 mjy, higher than for previous searches, probably because 
of the careful filtering of the candidates. 



Four of the 10 high-redshift QSOs were previ- 
ously known: 0747+27 (M. Irwin, private communication), 
0831+52 (Irwin et al 1998), 0839+51 (Snellen et al 2001) and 
1057+45 (Stern et al 2000). 0831+52 is unusually bright (E 
= 15) and the Lya forest been studied in detail (Ellison et 
al 1999). 327-MHz flux densities are available for three of 
the high-redshift QSOs from the WENSS survey (Rengelink 
et al 1997). The 0.3 - 1.4-GHz spectral indices a (S v oc v a ) 
are all flat: -0.1 for 0725+37, 0.0 for 0839+51, and -0.4 for 



1639+43. The other 7 QSOs are too faint to be detected in 
WENSS, or lie below the WENSS declination limit. 

The SEDs of the high-redshift radio QSOs reported here 
are as varied as those found in optical searches. At least 
0941+51 and 0831+52 show metal-line absorption redward 
of Lya, probably associated with saturated Lya lines in 
the forest. One, 1309+57, shows very narrow Lya emission, 
FWHM 20 A. 1639+43 has almost no Lya emission, which 
is unusual (only ~ 1/60 of SDSS QSOs), and it is similar in 
this respect to the z — 4.2 radio QSO, 0918+06, discovered 
by Snellen et al. Snellen et al suggested that the unusually 
strong absorption of the Lya line might be due to the host 
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galaxy (in both cases this is probably a giant elliptical). A 
few of the QSOs have probable damped Lya systems (DLAs, 
seen in ~ 20% of high-redshift QSOs). 0941+51 has a DLA 
at z = 3.52, and this has now been observed at high spectral 
resolution with WHT ISIS (Centurion et al, in preparation). 

The other 45 candidates are a mixed bag (see Maglioc- 
chetti et al 2000 for the typical distribution of identification 
types of all colours at 1 mjy). 1516+43 is a z — 2.59 broad- 
absorption- line (BAL) QSO with broad, deep SilV and CIV 
absorption troughs blueward of the emission features (Fig. 
5). 1035+30 (probably), 1114+29 and 1215+25 are low- 
redshift galaxies exhibiting narrow emission lines (Fig. 5). 
18 have SEDs consistent with those of giant elliptical galax- 
ies at z ~ 0.3. 23 have spectra of late-type M stars. This 
number is consistent with that expected for chance coinci- 
dences, but a few may be true radio stars (Helfand et al 1999 
describe a search for counterparts of bright stars in FIRST) . 



5 CONCLUSIONS 

We have conducted a pilot search for high-redshift radio 
QSOs, obtaining optical spectra of 55 FIRST radio sources 
(Si. iGHz > 1 mjy) with red (O-E > 3) starlike optical iden- 
tifications. 10 of the candidates are QSOs with redshifts 3.6 
< z < 4.4 (4 previously known). Six have z > 4. The re- 
maining 45 candidates are a mixture of low-redshift galaxies 
(misclassified as stellar) and M stars (mainly misidentifica- 
tions). The success rate (high-redshift QSOs / candidates) 
for this search is 1/2 for Siagh* > 10 mjy, and 1/9 for 
Siaghz > 1 mjy. With an effective search area of 4030 
deg 2 , the surface density of z > 4 QSOs discovered with 
this technique is 0.0015 deg -2 . 
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